Block Spin Ground State and 3-Dimensionality of (K,Tl)^Fei 6Se2 
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The magnetic properties and electronic structure of (K,Tl)yFei.6Se2 is studied using first-principles calcu- 
lations. The ground state is checkerboard antiferromagnetically coupled blocks of the minimal Fe4 squares, 
with a large block spin moment ^ 11.2/iB. The magnetic interactions could be modelled with a simple spin 
model involving both the inter- and intra-block, as well as the n.n. and n.n.n. couplings. The calculations 
also suggest a metallic ground state except for ^ = 0.8 where a band gap ^ 400 — 550 meV opens, show- 
ing an antiferromagnetic insulator ground state for (K,Tl)o.8Fei.6Se2. The electronic structure of the metallic 
(K,Tl)^Fei.6Se2 is highly 3-dimensional with unique Fermi surface structure and topology. These features 
indicate that the Fe-vacancy ordering is crucial to the physical properties of (K,Tl)yFe2-ccSe2. 

PACS numbers: 71.10.Hf, 71.27.+a, 71.55.-i, 75.20.Hr 



Superconductivity (SC) with moderate high transition tem- 
peratures [1-6] has been observed in a broad family of the 
iron-based materials. They are typically represented by the 
nil-type LaFeAsO yj, the 122-type BaFe2As2 El, the 
Ill-type LiFeAsli, and the 11-type FeSe ll. The par- 
ent compounds of these materials show a universal strip-like 
(collinear) antiferromagnetic (SDW) order ^lo', except for 
the 1 1 -type iron chalcogenides where the magnetic order is bi- 
collinearll^jisl]. The magnetic properties are closely related 
to a common two-dimensional Fe-atom square lattice and the 
electronic structures are featured by the cylinder-like hole and 
electron pockets around the F and M points respectively, with 
relatively weak dispersions alone the c-axis. By electron or 
hole doping the Fermi surfaces evolve smoothly in accordance 
with rigid band shift and the SC instability is enhanced once 
the magnetic order is suppressed [ 14]. 

Recently, a new family of the 122-type FeSe compounds 
K^Fe2Se2 iQ and (Tl,K)^Fe2-xSe2 flR H have been 
found to exhibit SC with transition temperatures Tc ~ 30 K. 
Moreover, the iron deficient compound (Tl,K)^Fe2-xSe2 
shows two remarkable features: (i) The SC (appears for x ~ 
0.12 — 0.3, y ~ 1) is in proximity to an insulating phase 
( for relatively larger x, or y ^ 0.8)1^-^8], and (ii) the 
Fe- vacancies may exhibit some ordered superstructures [Hill . 
Early Mossbauer experiment for TlFe2-xSe2 ill and re- 
cent transmission electron microscopy on KFe2-a^Se2 (for 
x = 0.4~0.5)|20] provide clear evidence for the tetragonal 
and orthorhombic superstructures in the FeSe layer. Previ- 
ous first-principles calculations suggested that the Fe-vacancy 
orthorhombic superstructure could be stablized with an stripe- 
like (collinear) AFM ground state in fTl.K^Fei ^Se9 i2lL l22h. 
The insulating behavior with an activation gap ~ 60 meV in 
transport observed for x ~ 0.5 or around[16] could be at- 
tributed to a moderate large short-ranged Fe-3d electron corre- 
lation, manifesting a possible Mott insulator driven by kinetic 
energy reduction due to the ordered Fe-vacancies[22]. The 
Mott-transition can be indeed realized by a relatively smaller 



Uc in a two-orbital model with vacancy orderings [l23i 12411 . 



So far it is yet unclear about the precise location of the crit- 
ical Fe-deficiency Xc (for y = 1) where the metal-insulator 
transition takes place. One should notice that in experiments 
the real Fe-content is sample dependent and may deviate from 
the nominal ones, and possibly ( in addition to introduce elec- 
trons to the FeSe layer), the role of Tl atoms is to stablize the 
Fe-vacancy orderings while the role of alkaline atoms (K, Rb, 
Cs) is to achieve higher Fe-content. Among all the iron defi- 
cient compounds, (K,Tl)2;Fei.6Se2 (x = 0.4) is of special in- 
terest due to its closer proximity to the transition point Xc and 
the peculiar ordering pattern of the Fe- vacancies which can 
be stablized in tetragonal crystalline superstructure 0161 [2Q|i . 
This is the simplest vacancy superstructure with the highest 
symmetry since all iron atoms are 3 -coordinated equivalent ly. 
Thus it is especially interesting to understand the electronic 
and magnetic structures of (K,Tl)^Fei.6Se2with the tetrago- 
nal Fe-vacancy superstructure. 

In this paper, we performed extensive study on TlFei.6Se2, 
using the first-principles simulations. We have also per- 
formed calculations on KFei.6Se2 for the magnetic structure 
as well as the density of states (DOS), which agrees well with 
TlFei.6Se2 results. We found that tuning y only leads to the 
rigid band shift for a specific Fe-deficiency x. Thus our results 
are valid for the mixture system (K,Tl)2yFei.6Se2. In particu- 
lar, we used the Vienna Ab-initio Simulation Package (VASP) 
11251 [2fl], which employs the plane- wave basis set and the pro- 
jected augmented wave (PAW) method | 27]. A body-centered 
orthorhombic primitive cell (FIG. |l(a)| l was used throughout 
the calculation unless otherwise specified. A 360 eV energy 
cut-off and a 4 x 4 x 4 F-centered k-grid were chosen to ensure 
the convergence of the total energy to 1 meV/cell. All the ge- 
ometry were optimized until the forces on each atom smaller 
than 0.01 eV/A and the total pressure smaller than 0.5 kB. For 
the DOS calculations, a much finer k-grid of 16x16x16 and 
the tetrahedra method were used. 

Similar to (K,Tl)Fei.5Se2, (K,Tl)^Fei.6Se2 has three dif- 
ferent stacking patterns. Here we focus on the in-plane mag- 
netic structure and consider the A A- stacking only. As in- 
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FIG. 1 : 1 1(a) [ The geometry and primitive cell from top-view. The shaded (blue) region indicates the fundamental block with the four Fe-atoms 
at the comers. The area encircled by the dashed line for the AFM2 configuration (the ground state) is 74.75 A^. The periodic boundary 
condition for such blocks extended over the whole lattice is impose d. 1 1(b) [ The proposed magnetic couplings. (Ji, J2) and (J{, J2) represent 
the intra-block and inter-block (n.n., n.n.n.) couplings, respectively. 1 1 (c) | to Various magnetic configurations. The red/blue atoms indicate 
the Fe-atoms with positive/negative total magnetic moment, respectively. For the AFMl configuration, two primitive cells consist a magnetic 
unit cell. In all figures, we show only the Fe atoms to enhance the visibility. 



dicated in (K,Tl)Fei.5Se2, the stacking ordering contributes 
only a negligible secondary correction to the total energy un- 
less it changes the symmetry of the crystal lattice! 22]. We also 
consider the anti-ferromagnetic (AFM) inter-layer coupling 
while its magnitude is negligible owing to the large inter-layer 
distances. As a check, we performed test calculations on one 
of our spin configurations. It turns out that inter-layer mag- 
netic coupling contribution is < 1 meV/Fe for TlFei.6Se2. 

To begin with, we first study the ground state magnetism 
by considering 8 possible in-plane AFM configurations [FIG. 
|l(c)| to p^)| , as well as the non-magnetic (NM) and ferro- 
magnetic (FM) orderings. We list the relaxed geometry 
parameters as well as their relative energies in TABLE H 
The AFMO configuration could be regarded as the checker- 
board AFM; whereas the AFMl and AFM6 configurations are 
the bi-collinear and zig-zag collinear orderings respectively. 
Our calculations suggest a ground state of the AFM2 type, 
whose configuration energy is 433 meV/cell (or 54 meV/Fe) 
lower than the second lowest (AFM6) configuration. Due to 
the symmetry of Fe- vacancies, all Fe sites are equivalent in 
(K,Tl)^Fei.6Se2, forming perfect square Fe blocks (indicated 



TABLE I: Lattice constants and magnetic properties of TlFei.6Se2. 
The lattice constants are transformed to represent 122 crystal. 





a{b) (A) 


c(A) 


rriFe (Ms) 


Ea (meV/Fe) 


NM 


3.8649(3.8649) 


13.1812 








FM 


3.7903(3.7904) 


14.7087 


2.8 


-62 


AFMO 


3.9026(3.9026) 


13.6949 


2.3 


-131 


AFMl 


3.8494(3.8182) 


14.1779 


2.7 


-183 


AFM2 


3.8667(3.8668) 


14.2420 


2.8 


-254 


AFM3 


3.8043(3.8517) 


14.1884 


2.7 


-199 


AFM4 


3.8892(3.8892) 


13.9058 


2.5 


-175 


AFM5 


3.8882(3.8778) 


13.9053 


2.6 


-180 


AFM6 


3.7645(3.8440) 


14.3045 


2.8 


-200 


AFM7 


3.8994(3.8385) 


14.0930 


2.6 


-183 



by the blue units in FIG. |l(b)| ) intercalated by Fe- vacancies. 
The AFM2 configuration can thus be regarded as checker- 
board antiferromagnetically coupled blocks of parallel aligned 
spins. Furthermore, due to the structure distortion induced by 
the Fe-vacancies, the magnetic couplings could be classified 
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into two groups: the ones within a square block (intra-block) 
and the ones across two near-by square blocks (inter-block) . 
To model our energetic results, we incorporated a spin model 
involving both the nearest neighbour (n.n.) and the next near- 
est neighbour (n.n.n.) couplings: 

n,(y. 

Where, n denotes the block index, 5 is short for the nearest 
neighbouring block to block n, a is the site-index which goes 
from I to 4, as selects the site connecting to the nearest neigh- 
bouring block J; Ji and J{ ( J2 and J2) indicate the n.n. (the 
n.n.n.) couplings of intra- and inter-block, respectively. If we 
further use the approximation that only the 5*^ component is 
involved (Ising model with being the same in all configu- 
rations at each Fe sites), we could fit the energetics of the 8 
AFM configurations using the least squares method to obtain 
J = 2JS'^. The resulting intra-block couplings Ji and J2 
are -86 meV and -9 meV, respectively; while the inter-block 
couplings J[ and J2 are -29 meV and 38 meV, respectively, 
with fitting correlation a = 97.04%. The AFM inter-block 
n.n.n. interaction J2 dominates the inter-block interactions 



and both intra-block interactions are FM, thus a block-type 
checkerboard AFM configuration is favored. 

We then examine the electronic structure of 
(K,Tl)^Fei.6Se2. We show the band structures of the 
AFM2 and NM states for ^=1.0 in FIG. [2(b)] and [2(i)l 
respectively. The band structure of the AFM2 state suggests a 
metallic nature, with a band gap ~ 400 meV for TlFei.6Se2 
(or 550 meV for KFei.6Se2) for the states 16 meV (or 39 
meV ) below Ep. Interestingly, the DOS from the top of 
the band gap to Ep integrates to exactly 0.2 electron per 
(K,Tl)^Fei.6Se2 (^=1.0) formula, suggesting that the material 
would become an insulator if (K,T1) content is decreased 
by 20%. As a check, we have further performed calcula- 
tions for K0.8Fe1.6Se2, assuming two types of K- vacancy 
orderinss 1I32I1. Then a band gap ~ 600 meV shows up at 
for either cases. The result implies that the K- vacancies shift 
the chemical potential but their orderings do not change the 
band structures. 

Another interesting feature of the (K,Tl)2yFei.6Se2 elec- 
tronic structure is that in the NM state the band energies dis- 
perse significantly along /c^-axis, manifesting its highly 3- 
dimensional characteristics. This feature is fundamentally dif- 
ferent from the KFe2Se2 compound[28] and all other iron 
pnictides. Of course, both the NM and AFM DOS indicate 
that the Fe-3d orbitals dominate the states near Ep, similar to 
all iron-based superconductors. 



4 











Ti/a 



(c) 



(d) 



FIG. 3: Femi surfaces of TlFei.6Se2reconstructed using MLWFs 
at (a,c): non-magnetic (NM) state and (b,d): block- spin anti- 
ferromagnetic state; (c,d): Cross section of the Fermi surface of 
TlFe2-xSe2at kz = 0.0 plane with x = 0.4. Both NM and block 
spin AFM Fermi surfaces are 3D-like. 



Finally, we reconstruct the Fermi surfaces (FIG. |3} of 
the metallic (K,Tl)^Fei.6Se2 (^=1.0) by fitting the band 
structure using the maximally localized wannier functions 
(MLWFs) 111, til. The Fermi surface of both NM and block- 
spin AFM (K,Tl)^Fei.6Se2 are highly 3- dimensional, al- 
though their specific shape is quite different from each other. 
It is worthwhile to notice that due to the Fe- vacancy super- 
structure the first Brillouin zone (BZ) is not the same as 
the one in (K,Tl)Fe2Se2, and thus the M point in FIG. [5] 
is not (7r/a,7r/a). Nevertheless, the Fermi surface topology 
of the present (K,Tl)^Fei.6Se2 compound is quite unique 
compared to either KFe9Se9 ll28ll or all other iron pnictides. 
It strongly indicates that the formation of Fe-vacancy su- 
perstructures is crucial to the electronic structures of the 
(K,Tl)^Fe2-a;Se2 compounds. Unlike the change in y, the 
topological change of the Fermi surface across x ^ OA indi- 
cates that the electronic and magnetic structures for different 
Fe-vacancy ordered materials can not be approached by a rigid 
shift of the chemical potential. 

It is worth noting that while the Fe-vacancy ordering is 
crucial to the block spin magnetic pattern and the finite gap 
for (K,Tl)o.8Fei.6Se2, the strong disorder of the Fe- vacancies 
may destroy the magnetic ordering leading to the metallic 
ground state. In the vacancy disordered state, the first Bril- 
louin Zone would remain the same as (K,Tl)Fe2Se2 and thus 
the rigid-band model would remain effective. The whole elec- 
tronic structure in that case could then be approximated with 
a hole-doped (K,Tl)^Fe2Se2l28]. Actually, it has been sug- 
gested that randomly distributed Fe- vacancies in the doped an- 



tiferromagnetic (or Mott) insulator may lead to a spin- sing let 
s- wave superconductor [24J . 

In conclusion, we have performed first-principles calcula- 
tions on (K,Tl)2yFei.6Se2. A block-type checkerboard anti- 
ferromagnetic ground state was identified and the AFM inter- 
block n.n.n. coupling interaction dominates. Our calculations 
suggest a metallic ground state for y = 1 with a 400-550 meV 
band gap which appears slightly below the Fermi level and 
an insulating ground state for (K,Tl)o.8Fei.6Se2. The exper- 
imentally observed insulating behavior may be due to both 
the 20% (K,T1) deficiency and the iron vacancy superstruc- 
ture. Furthermore, the electronic structures of the metallic 
states show a significant 3 -dimensional feature with a unique 
Fermi surface topology, indicating that the formation of Fe- 
vacancy superstructure is crucial to the physical properties of 
(K,Tl)^Fe2-^Se2 . 
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Note added: Recently, we become aware of a paper by 
W. Bao et al [31] on neutron diffraction experiment for 
K0.8Fe1.6Se2. Our result is in agreement with the reported 
magnetic ordering pattern. 
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